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The electromagnetic interactions of nuclei have been used extensively to elucidate their structure 
and dynamics. In the early days of nuclear physics, the magnetic moments of the light nuclei 
helped to reveal that they behaved like a collection of weakly interacting nucleons that, to a very 
large degree, retained their identity, despite being bound together by the strong nuclear force. This 
feature, in part, led to the establishment of the nuclear shell model as a phenomenological tool 
with which to predict basic properties of nuclei throughout the periodic table. The strong nuclear 
force emerges from quantum chromodynamics (QCD) as a by-product of confinement and chiral 
symmetry breaking. The fact that, at the physical values of the quark masses, nuclei are not simply 
collections of quarks and gluons, defined by their global quantum numbers, but have the structure 
of interacting protons and neutrons, remains to be understood at a deep level. In this presentation, 
I will show the results of recent Lattice QCD (LQCD) calculations of the magnetic moments [1] 
and polarizabilities [2] of light nuclei at pion masses of tuti; ^ 450 MeV and 805 MeV, along with 
the first calculation of the radiative-capture cross section of np dy [3] extrapolated to nature and 
compared with experiment. 

Lattice QCD is moving rapidly into a position to calculate the properties and interactions 
of the lightest nuclei at the physical values of the light-quark masses and with the inclusion of 
fully dynamical QED. The ultimate reason for such calculations is to refine the nuclear forces, 
and enhance our ability to reliably predict the properties and interactions of nuclei throughout the 
Periodic Table, beyond what is possible with experiment alone. Currently, we are only now en¬ 
tering the “verification stage”, where experimentally measured nuclear physics observables that 
can be accessed by LQCD are demonstrated to be reproduced within the uncertainties of the mea¬ 
surements and of the calculations. It will likely be another few years before calculations are of 
sufficient precision and accuracy to have verified LQCD as a reliable predictive tool for nuclear 
physics. Until recently, it was only the scattering of two nucleons [4, 5, 6, 7] and the binding 
energies of the lightest nuclei and hypernuclei [8, 9, 10, 11, 12, 13] that were being pursued with 
LQCD because calculations of other properties, such as magnetic and quadruple moments, are not 
meaningful without clearly establishing a bound nucleus. Multi-meson systems have been consid¬ 
ered previously as somewhat of a prelude to multinuclear systems, in particular the extraction of 
multi-body interactions [14, 15, 16, 17]. 

Conceptually, Lattice QCD calculations are a straightforward and brute-force evaluation of 
the QCD path integral for the observables of interest. However, this is not quite QCD that is being 
calculated as the evaluation is performed in Euclidean space with a finite lattice spacing providing 
the ultra-violet regulator and in a finite volume which modifies the infrared. Eor a lattice spacing 
that is small compare with typical strong interaction length scales, the Symanzik action is used to 
systematically remove lattice spacing artifacts, usually through performing multiple calculations 
of the same quantities over a range of lattice spacings. Eor a lattice volume that is large compared 
with the size of a compact state, such as a nucleus, effective field theory (EET) techniques can be 
used to extrapolate to infinite volume, by performing calculations in multiple volumes. Scattering 
states are a different matter, and the geometry of the finite-volume, along with the momentum of the 
system in the volume and the chosen boundary conditions, determines the nature of the eigenstates 
and their relation to Minkowski-space S-matrix elements. One should view LQCD as technique 
to solve a low-energy EET of quarks and gluons that is applicable for momenta that are much 
below the inverse lattice spacing of the lattice. Typical lattice spacing used in present-day LQCD 
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calculations are 0.1 fm, corresponding to a momentum scale of ^ 2 GeV. 

One of the important features of LQCD calculations of systems involving baryons (compared 
to mesons) is the signal-to-noise degradation that is encountered in the correlation functions. At 
large times, the signal-to-noise ratio in a nucleon correlation function scales as 
degrades somewhat less severely at intermediate times by judicious choices of interpolating opera¬ 
tors. Insight can be gained into the behavior of correlation functions from the results of calculations, 
such as those shown in Fig. 1. The fluctuations in the gauge field clearly impact the evolution of 



Figure 1: The same spatial slice of the correlation functions of the n, A and AA at three different times 
superposed on the action density of the gluon field. These “toy” calculations are quenched with a DBW2 
gauge action and a pion mass of mTu ^ 350 MeV. 

the hadron correlation functions in highly non-trivial ways, as one expects. These interactions in¬ 
troduce significant non-Gaussianity into the distribution of correlations functions on any particular 
time-slice. 

NPLQCD has performed a comprehensive study of s-shell nuclei and hypernuclei, and of 
nucleon-nucleon scattering, at the flavor-SU(3) synunetric point at pion mass of ruji ^ 805 MeV, 
where the strange-quark mass is tuned to its value in nature. Calculations were performed in 
multiple lattice volumes, parametrically reducing the residual finite-volume effects. However, due 
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to computational resource limitations, calculations were performed at one lattice spacing only, but 
one that is estimated to be small enough to have only small effects on the values of calculated 
binding energies and scattering parameters. The calculated binding energies are shown in Fig. 2. 
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Figure 2: The spectrum of the s-shell nuclei and hypernuclei obtained by the NPLQCD collaboration [11] 
at the flavor-SU(3) symmetric point with a pion mass of ^ 805 MeV. 



Ground states, and in some cases excited states, are well resolved for nuclei and hypernuclei with 
atomic number A = 1,2,3 and 4. Isolating the ground states is necessary for further studies of 
their properties, such as magnetic moments and polarizabilities. The scattering phase shift in the 
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Figure 3: Nucleon-nucleon scattering phase shifts (left panel) and ^cot5 (right panel) in the -channel at 
the flavor-SU(3) symmetric point with a pion mass of ^ 805 MeV [7]. 

^S\ channel, shown in Fig. 3, is extracted from the location of higher-lying states in the spectrum 
using Luscher’s quantization conditions (and assuming higher partial-wave phase shifts are small). 
Also shown in Fig. 3 is the extracted values of ^cot5 that uniquely defines the nucleon-nucleon 
scattering amplitude below the start of the t-channel cut. The precision of the calculation allows 
for a determination of both the scattering length and effective range. The scattering length is 
found to be twice as large as the effective range, suggesting that the deuteron remains a somewhat 
large object over a significant range of quark masses, and is not finely-tuned as had been thought 
previously. This intriguing feature requires further investigation. 
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When the results from NPLQCD [9, 11,18] are combined with those from Yamazaki et al. [12, 
13], a pattern begins to emerge as to how the deuteron binding energy varies with the light-quark 
masses, as shown in Fig. 4 ^ It appears that nothing dramatic is going to happen to the deuteron 
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Figure 4: Summary of fully-dynamical Lattice QCD calculations of the deuteron binding energy as a func¬ 
tion of the pion mass [9, 11, 12, 13, 18]. 


binding as the light-quark masses are reduced to their physical values, and that the system evolves 
smoothly to the experimentally measured binding energy. Calculations at a pion mass of niji ^ 
300 MeV are ongoing, and we expect that calculations closer to niji ^ 140 MeV will be performed 
soon, see Ref. [22]. 

Having identified ground state plateaus of the lightest nuclei, we are now in a position to 
ask questions about their structure. Historically, the magnetic moments of nuclei, and magnetic 
transitions between states, have revealed much about their structure. In particular, that they are 
close to the values predicted by the naive nuclear shell model. Exploring the range of validity of 
the nuclear shell model as a function of the light-quark masses - a dimension that is not accessi¬ 
ble experimentally - may reveal further secrets about the nature of the strong interactions and the 
structure of nuclei. In fact, it does. There are a couple of techniques with which to determine the 
magnetic moments of hadrons with LQCD calculations, and NPLQCD chose to use background 
magnetic fields in its studies. A C/(l) gauge field is constructed to produce a uniform and time- 
independent magnetic field (along the z-direction for simplicity). It is required to be quantized (in 
units of h), as shown in Fig. 5. Each configuration of QCD gauge fields in a given ensemble is 
post multiplied by this U{1) gauge field, which subjects the valence quarks to the magnetic field, 
but not the sea quarks. Having the sea quarks also interact with the magnetic field would require it 
to be present during the generation of the ensemble, which is currently beyond our computational 

Tt should be mentioned that the calculations of the HAL QCD collaboration do not find bound nuclei at these 
heavier pion masses, e.g. Ref. [19]. However, they have not directly calculate bindings, but have used a uncontrolled 
truncation of a non-local lattice correlation function to arrive at their conclusion(s). The apparent disagreement is likely 
due to the limitations of this hadronic modeling. For more discussion on this technique, see Refs. [20, 21] . 
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Figure 5: A f/(l) gauge field that generates a uniform and time-independent magnetic field 


resources. Therefore, the results of such calculations are generally not predictions of QCD, due 
to the absence of electromagnetic coupling to the sea quarks, however there are certain quantities 
that are predictions of QCD. For instance, at the flavor-SU(3) symmetric point, the electric charge 
matrix is traceless, and as such there is no coupling of one insertion of the magnetic field to the 
sea quarks. Therefore, the magnetic moments of the hadrons and nuclei calculated with such a 
post-multiplication by the background field are equivalent to those generated with fully-dynanfical 
QED up to terms that are higher order in In contrast, while the isovector polarizabilities, cor¬ 
responding to two insertions of the magnetic field, are equivalent to a full calculation, the isoscalar 
polarizabilities are not due to loop diagrams involving two photons coupling to one sea-quark loop. 
The complexity is more significant away from the SU(3) point as loops involving one photon no 
longer vanish by the tracelessness of the charge matrix due to SU(3) breaking effects. 

A complication of LQCD calculations in the presence of a background magnetic field is that 
the energy eigenstates of charged states are no longer momentum eigenstates, but are Landau levels 
in the plane transverse to the field and momentum eigenstates along the direction of the field. 
In dimensionless units, the energy difference of the ground state of a nucleus, h, of spin j and 
projection along a magnetic field along the z-direction (B = and P|| = 0) is 


2^ 


aSEh-j^ = \ a^Ml + {2nL + l)Qha^e B\-aMh - —flhjza^\eB\ 
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where ni is the integer associated with the Landau level, jlh = is the dimensionless magnetic 
moment, both and 5/, are terms higher order in the expansion (but which are present to stabilize 
fits to the magnetic moments and polarizabilities). The dimensionless polarizability and tensor 
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polarizability are 

R _^N -Mn) (MO) ^ (2) _ Ml (Ma -Mn) ^ (Ml) 

Ph — ^2 Ph ^ Ph ~ g2 Ph ’ 

and a^\e B\ — ^ is the dimensionless field strength. 

In the absence of a magnetic field, nuclear correlation functions are constructed from “hadronic 
blocks”. These are three light-quark propagators contracted at the sink with the quantum numbers 
of a proton or neutron, and then momentum projected. This method of constructing nuclear corre¬ 
lation functions has proven to be efficient and successful for light nuclei. However, in the presence 
of a magnetic field, the Landau levels associated with a proton are different from a nucleus, say the 
triton, as both the charge and the mass determine the orbit. Therefore, it was anticipated that the 
hadronic blocks might not be optimal for the production of nuclear correlaton functions in a mag¬ 
netic field. This was observed numerically, as shown in the ratio Z{B) /Z(0) for the triton ground 
state in Fig. 6. The calculations were most precise for the neutron and di-neutron systems for which 



n 

Figure 6: The ratio of overlaps of the hadronic blocks onto the triton ground state as a function of magnetic 
field [2]. 

there are no Landau levels. 

For the simplest case of the neutron, the negatively-shifted spin-down state depends essentially 
linearly on the magnetic-field strength out to very high magnetic fields, as shown in Fig. 7. For the 
fields used in this work, |^B| ^ 0.05/z GeV^, and therefore h—\ corresponds to |B| 10'*^ Gauss. 

In physical units, the energy shift in the h— 12 magnetic field is ISE ^ 400 MeV. Therefore, we 
see the lower-lying state depend linearly upon the magnetic field out to about r\j 10^0 Gauss. The 
upper level, on the other hand, shows extensive curvature setting in at quite low fields (on this 
scale), and has similarities with what one observes in avoided crossing in two-state systems (with 
the other state not observed in this case). It is clear from this behavior that the magnetic polariz¬ 
ability of the neutron is determined entirely by the curvature of the upper state, and receives little 
or no contribution from the lower state. An interesting question to ask, that obviously requires 
significantly finer lattices, is to how the states behave in the limit of asymptotically large magnetic 
fields, where one expects only the behavior of nearly massless charged quarks to dominate the 
spectrum. In particular, when does the expected curvature set in for the lower state? The behavior 
of the proton states is more complex due to the Landau-level structure. However, the system can 
be analyzed and the magnetic moment and polarizability extracted in a similar way to that of the 
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Figure 7: The energy shifts induced in each spin state of the neutron (left panel) and proton (right panel) by 
the magnetic field [2]. 

neutron. The uncertainties are larger as the correlation functions are of somewhat poorer quality 
due to smaller overlap factors. Further, the landau level associated with the plateau on the cor¬ 
relation function is unknown and must be deternfined by fitting [2]. A summary of the magnetic 



Figure 8: The results of the Lattice QCD calculations of nuclear magnetic moments [1] are shown as the 
blue bands, while the corresponding experimental values are shown by the red dashed lines. Natural nuclear 
magnets have been used, in which the nuclear magneton is defined in terms of the nucleon mass at the given 
value of the light-quark masses. 

moments of the nucleon and lightest nuclei is shown in Fig. 8. It is remarkable that, when given 
in terms of natural Nuclear Magnetons (defined with the nucleon mass at the given light-quark 
masses), the magnetic moments of the light nuclei at a pion mass of tuti; ^ 805 MeV are very close 
to their values at the physical light-quark masses. This implies that essentially all of the light-quark 
mass dependence of the magnetic moments is determined by the mass of the nucleon, and that a 
non-relativistic quark model type scenario is providing the dominant contribution - one in which 
a naive weighted sum of quark-model quark spins is dominant. This could entirely be a conse¬ 
quence of the large-Ac limit of QCD or it might be something more. Another, and perhaps more 
remarkable, feature is that the magnetic moments of the nuclei are essentially given by the sum 
of the nucleon magnetic moments in a naive nuclear shell-model configuration. For these nuclei, 
this is somewhat trivial compared to larger nuclei, but nonetheless the two neutrons in the triton 
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are largely in a spin-zero configuration with the magnetic momentum being essentially that of the 
proton. The deviations observed in nature from the naive shell-model prediction are in agreement, 
within uncertainties, with the result of the LQCD calculation at the heavier pion mass. This leads 
one to observe that nuclei behave as collections of weakly interacting nucleons over a large range 
of light-quark masses, and that the phenomenological nuclear shell-model is not limited in appli¬ 
cability to the physical point, but is somewhat of a generic feature of QCD for arbitrary light-quark 
masses. It will be interesting to learn if there are values of the quark masses for which nuclei col¬ 
lapse into a strongly interacting configuration of quarks and gluons rather than of weakly nucleons 
with a hierarchy of multi-nucleon forces. The magnetic moments have also been calculated at a 
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Figure 9: The isovector anomalous magnetic moment of the nucleon [from the NPLQCD collaboration]. 


pion mass of ^ 450 MeV. The same behavior is observed, as demonstrated in the isovector 
anomalous magnetic moment of the nucleon shown in Fig. 9. 

A key reaction in Big Bang nucleosynthesis is the radiative capture process np dy, which is 
the primary mechanism to first produce deuterium. It has been long known that a naive calculation 
of the cross section for this process using only nucleon-nucleon scattering S-matrix elements pro¬ 
vides only ^ 90% of the experimentally measured cross section at low energies. It has also been 
long recognized that this defficiency is due to the relatively large role played by meson-exchange 
currents (MFCs), or alternatively from a modern perspective the role of correlated short-distance 
two-nucleon interactions with the electromagnetic field. Once the contributions from these inter¬ 
actions are included, the cross section comes into agreement with experiment. ^ With the same set 
of LQCD calculations that provided the nuclear magnetic moments (at both = 450 MeV and 
805 MeV), the splitting between the = 0 neutron-proton energy levels in the background mag¬ 
netic field allowed for a determination of the correlated short-distance two-nucleon interactions 
with the magnetic field, and hence for the first QCD prediction of np dy [3]. At the heavier pion 
mass, this enabled a prediction of the cross section without any experimental input - truly a QCD 
prediction, while at the physical pion mass, the scattering parameters were used, in part, to predict 

^In low-energy EFTs, there is an Ml local operator with an coefficient that cannot be determined from scattering 
alone [23, 24], and hence agreement at one kinematic point is somewhat trivial. However, after precisely determining 
this counter term, there is good agreement in the Ml amplitude over a range of energies [25, 26, 27] that would otherwise 
have been in confiict. 


9 





Magnetic properties of light nuclei from lattice QCD 


Martin J. Savage 


the complete cross section. Using the physical scattering parameters and the LQCD calculations 
of the correlated short-distance two-nucleon interactions with the magnetic field, a cross section 
of = 332.4^4'^7 mb is calculated at a neutron incident speed of v = 2200 m/s, which is to be 
compared with the experimental value of = 334.2 zb 0.5 mb. 

The curvature of the energy of the nucleon or nucleus, after removing the contribution from 
the Landau level, provides a determination of its magnetic polarizability. In nature, there is signif- 
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Figure 10: The magnetic polarizabilities of the lightest nuclei at a pion mass of m^j; ^ 805 MeV [2]. The 
left panel is in units of the naive A-pole contribution, e^/M'^{M a—M^), while the right panel is in physical 
units. 


icant cancellation between the contribution from the A-pole and from chiral loops to the magnetic 
polarizability, a cancellation that is expected to diminish as the quark masses are increased [28] 
(the A-nucleon mass splitting is relatively insensitive to the light-quark mass). A large isovector 
component to the nucleon magnetic polarizability is found, and as mentioned previously, while 
the isoscalar polarizabilities are subject to modifications (that are expected to be small) due to the 
omission of disconnected diagrams, the isovector contributions are complete at the flavor-SU(3) 
symmetry point. It is interesting that, as shown in Fig. 10, the magnetic polarizabilities of the light 
nuclei are found to be near that of the proton. 

The precision with which we have been able to determine the neutron systems is sufficient to 
determine that, while the di-neutron is bound at these heavier quark masses, there are values of the 
magnetic field for which it unbinds and the ground state becomes two isolated neutrons. This is a 
QCD Feshbach resonance! Figure 11 shows the results of the LQCD calculation at = 805 MeV 
of the energy difference between the bound di-neutron and two isolated neutrons. Clearly, the 
di-neutron is becoming less bound with increasing magnetic field, and is consistent with being 
unbound, and hence the two-neutron system having an infinite scattering length, near h ^ 5. 

In summary. Lattice QCD is emerging from an extended research and development phase into 
the production phase. Calculations of the binding and properties of light nuclei are now possible, 
and I have presented the state-of-the-art of such calculations. The magnetic moments are providing 
important insights into the nature of nuclei and their stability with regard to changes in the fun¬ 
damental parameters of nature. The first inelastic nuclear reaction, np dy, has been calculated 
and when extrapolated to the physical point is found to be in agreement with experiment. The 
magnetic polarizabilities of light nuclei have been calculated, and a large isovector component has 
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Figure 11: The energy difference (in lattice units) between two neutrons and the bound di-neutron as a 
function of magnetic field strength [2]. The red band corresponds to breakup threshold. 

been identified in the nucleon. These works are merely a peek into the future precision era of QCD 
calculations of low-energy nuclear physics observables. 

Calculations were performed using computational resources provided by the Extreme Sci¬ 
ence and Engineering Discovery Environment (XSEDE), which is supported by National Science 
Eoundation grant number OCI-1053575, NERSC (supported by U.S. Department of Energy Grant 
Number DE-AC02-05CH11231), and by the USQCD collaboration. This research used resources 
of the Oak Ridge Leadership Computing Eacility at the Oak Ridge National Laboratory, which 
is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE- 
AC05-00OR22725. The PEACE Research Infrastructure resources Curie based in Erance at the 
Tr'es Grand Centre de Calcul and MareNostrum-111 based in Spain at the Barcelona Supercom¬ 
puting Center were also used. Parts of the calculations used the Chroma software suite [29]. This 
work was supported in part by DOE grant No. DE-EG02-00ER41132. 
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